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ABSTRACT 



Context. Several studies have tried to ascertain whether or not the increase in abundance of the early-type galaxies (E-SOa's) with 
time is mainly due to major mergers, reaching opposite conclusions. 

Aims. We have tested it directly through semi-analytical modelling, by studying how the massive early-type galaxies with 
log(Al,/A4o) > 11 at z ~ (mETGs) would have evolved backwards-in-time, under the hypothesis that each major merger gives 
place to an early-type galaxy. 

Methods. The study was carried out just considering the major mergers strictly reported by observations at each redshift, and assum- 
ing that gas-rich major mergers experience transitory phases as dust-reddened, star-forming galaxies (DSFs). 

Results. The model is able to reproduce the observed evolution of the galaxy LFs at z < 1, simultaneously for different rest-frame 
bands (S, /, and K) and for different selection criteria on color and morphology. It also provides a framework in which apparently- 
contradictory results on the recent evolution of the luminosity function (LF) of massive, red galaxies can be reconciled, just consider- 
ing that observational samples of red galaxies can be significantly contaminated by DSFs. The model proves that it is feasible to build 
up ~ 50-60% of the present-day number density of mETGs at z < 1 through the coordinated action of wet, mixed, and dry major 
mergers, fulfilling global trends that are in general agreement with mass-downsizing. The bulk of this assembly takes place during 
~ 1 Gyr elapsed at 0.8 < z < 1, providing a straightforward explanation to the observational fact that redshift z ~ 0.8 is a transition 
epoch in the formation of mETGs. The gas-rich progenitors of these recently-assembled mETGs reproduce naturally the observational 
excess by a factor of ~ 4-5 of late-type galaxies at 0.8 < z < 1, as compared to pure luminosity evolution (PLE) models. 
Conclusions. The model suggests that major mergers have been the main driver for the observational migration of mass from the 
massive-end of the blue galaxy cloud to that of the red sequence in the last ~ 8 Gyr. 



Key words, galaxies: elliptical and lenticular, cD — galaxies: evolution - 
luminosity function, mass function 
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1. Introduction 

. I_| ' Early-type galaxies (ETGs) have become cornerstones in our un- 
. derstanding of the mass assembly in the Universe, as they domi- 
j_i ■ nate the massive end of the galaxy mass function at z < 0.8, host- 
ed ' ing more than a half of t he total local stellar content (Bel l et alj 
' I2004bl lllbert et alj|20l0h . Although hierarchical models are the 
most successful in reproduc ing the general properties of the 
Universe at the present (see iNavarro & Steinmetzl [2000), sev- 
eral observational result s related to the formation of the mETGs 
call it into question (ICowie et al. _ 1996t 



2.5, see Eggen et al 



Kauffmann & White 



1962 



1993 



Heavens et all 12004: 



Bund v et alj2006tlJimenez et al .1120071: iDavies et al.l l2009). 



In particular, hierarchical and primordial-collapse theories of 
galaxy formation derive extremely-different assembly epochs for 
the mETGs. Actual ACDM models predict that these systems 
are the final remnants of the richest merging sequences in the 
Universe, and thus, the latest to be completely in place into the 
cosmic scenario (at z < 0.5), whereas monolithic collapse the- 
ories push their formation towards much earlier epochs (at z > 



; Too mre1[l977l . IWhite & Reesl [19781: 
Springel et al.l 12005b . In this sense, 
the observational phenomenon of galaxy mass-downsizing has 
become a challenge for hierarchical theories. Consisting on 
that the most massive galaxies (A1» > 5 x 1 n Alo) seem 
to ha ve have been in p lace since z ~ 2 - 3 dKodama et al.l 
120071: iDrorv & Alvarez! 120081: iKang & Iml 120091) , whereas the 
less-massi ve systems get their ac t ual volume densities at later 
epochs dGlazebrook et alJ l2004t iPerez-Gonzalez et al.l 2005, 
2008; Walche7etani008), mass-downsizing seems to favour 
a monolithic collapse origin for mETGs. This scenario is also 
supported by the substantial amounts of m ETGs that are be 
ing detected at z > 2 in the last years (e.g ., lElston etaUll988 
iFrancis et al. 



ICimatti et al 



\l99$L iHughes et al.lfT99l IWaddington et al.i 2002 
2004l: lFontanaetal.ll2009l) . and bv the negligible 



number evolution that several authors report for them up to 
z ~ 1 .2 dPozzetti et all 120031; iNolanl 120041; ICimatti et alj|2006t 
ISaracco et al.ll2009b iFerreras et alJl2009albh . 
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However, there are conflicting views on the amount of this 
number evolution, as the fraction of mETGs assembled since 
z ~ 1 ca n vary between ~ 20% and ~ 60%, depending on 
the author (Bell et al. 2004b; Ilbert et al. 2006; Faber et al. 2007; 



Scarl ata et alj I2007ac iCool et alj I2008t iMatsuoka & Kawaral 
2oToT F07 hereandafter). Moreover, the traces of past inter- 



actions exhibited by most of these systems and the merger 
fractions derived up to z ~ 1 suggest that major mergers 
must have contribute d significantly to the evolution of mETGs 
in the last ~ 8Gyr 



Conselice et al 



Conselice et al 



2003 



dLe Fevre et all 120001: iPatton et al. 



Ivan Dokkum 2005; Feldmann et al 



2002; 



2008 



20091) . This is also supported by recent observa- 
tions that prove the existence of a link between the AGN activity 
in low-re dshift ETGs and a significant merger event in their re- 
ce nt past dSchawinski et"aT]|2007L 1201 Oh . 

Ivan Dokkum et alJ (120001) showed that ETGs at high-redshift 
seem to form a homogeneous old population just because the 
progenitors of the youngest present-day ETGs are usually not 
included in the samples under study. This effect, known as pro- 
genitor bias, can extremely affect the conclusions derived on the 
evolution of the ETGs in those studies that exclude late-type 
progenitors f r om th eir samples. In fact, a pioneering study by 
iKavirai et alj ((2009a) has proven that, when the true progenitor 
set of the present-day ETGs is considered in the standard ACDM 
framework (regardless of their morphologies), less than 50% of 
the stellar mass which ends up in ETGs today is actually in early- 
type progenitors at z ~ 1, in agreement with recent observations. 
This result means that observations have not established yet how 
or even when the mETGs formed, nor if their apparently un- 
changed properties until z ~ 1 are really due to their evolution 
or to observational uncertainties and progenitor bias. 

Several studies have tried to ascertain whether or not the in- 
crease in abundance of ETGs with time arises primarily from 
major mergers, deriving opposite conclusions dBlantonl 2006; 
Bund v et all [20071 iDrorv & Alvarez! l2008t iBundv etalJ 12009: 
Cattaneo et alJl2010h . Most of these studies use merger trees 



from ACDM cosmological simulations, and some of them prove 
the feasibility of reproducing obs ervations within a theoret- 
ical ACDM framework (see, e.g., iKhochfar & Burkeril [2003; 
lLanzoni et"aU 2005; Ka virai et alj |2009a). But the real contribu- 
tion of the observed major mergers to the evolution of the high- 
mass end of the galaxy LFs is still unknown. So, we propose 
to test the hierarchical origin of mETGs by studying how the 
population of mETGs would have evolved backwards in time, 
under the hypothesis that they derive from the major mergers 
that are strictly reported by observations at each redshift. We 
intend to test if major mergers are really the main contributors 
to the buildup of mETGs or if additional processes are required 
to explain it. This is the first paper of a series summarizing the 
main results of this model. In the present paper (paper I), we 
describe the model and compare its predicted LFs with obser- 
vational data. Paper II analyses in detail the role of the obser- 
vational major mergers (wet, mixed, and dry) in the buildup of 
mETGs since z ~ 1 . 

The present paper is organized as follows. Section $2] is 
devoted to the model description. In $3] we compare the evo- 
lution of the LFs predicted by the model for different galaxy 
types with real data, using different galaxy-selection criteria (on 
color, spectral type, morphology). The compatibility of the re- 
sults with the mass-downsizing scenario and with other recent 
observational results is discussed in $4] In this section, we also 
comment on the model limitations. Finally, a brief summary of 
results and some conclusions are addressed in ^5] We will use 
a Q.m = 0.3, £2a = 0.7, Ho = 70 km s _1 Mpc~' concordant 



cosmology thr oughout the paper dNavarro &Steinmetzl l2000: 
IXu et al.l l20071. All magnitudes are given in the Vega system. 



2. The model 

The mod el adopts the b ackwards-in-time technique first intro- 
duced by Tinsley (1980). It traces back-in-time the evolution of 
the local galaxy populations considering two different sources 
of evolution: 1) the number evolution that observational merger 
rates imply at each redshift, and 2) the typical luminosity evo- 
lution (L-evolution) of each galaxy type due to its star forma- 
tion history (SFH). The evolution of the volume element derived 
from the assumed cosmology is also considered. All model pa- 
rameters and assumptions are strictly based on robust observa- 
tional and computational results. The main novelty of the model 
is the realistic treatment of the effects of major mergers on the 
LFs. It considers the different phases and time-scales in a major 
merger, as stated by observations and N-body simulations. As 
the fate of a merger remnant depends strongly on the gas-content 
of the progenitor galaxies, we have also had into account the dif- 
ferent properties of these progenitors depending on whether the 
merger is dry (between gas-poor galaxies), wet (both progeni- 
tors are gas-rich), or mixed (between a gas-poor galaxy and a 
gas-rich one), as well as the relative fraction of each merger type 
at each redshift according to observations. Two key assumptions 
relative to major merging are adopted in the model: 1) the fact 
that a major merger generates an E-SOa galaxy, and 2) that the 
gas-rich galaxies involved in a major merger experiment transi- 
tory phases in which they become DSFs. 

We use observationally-determined merger fractions to ac- 
count for the number density of merger events at each redshift. 
The merger fraction at a certain redshift is defined as the number 
of mergers taking place per galaxy at that redshift. Considering 
the volume enclosed at the considered redshift interval, we can 
estimate the number density of major mergers at redshift z, 
m (z), as: 

0m(z) = fm(z) ■ 0total(z), (1) 

being f m (z) the merger fraction and t otai(z) the total number den- 
sity of galaxies at z (i.e., the total galaxy LF at that redshift). 

2.1. The original ncmv code 



We have used the code NCMOD, created by iGardnerl (I1998L 
G98 hereandafter) for generating predictions on galaxy num- 
ber counts, as a basis for the model. The author kindly made 
the code available at his web site, and permited its use, 
modification, and distribution without charge for the purposes 
of scientific research. Briefly, the code evolves the local LF 
back-in-time, for a number of galaxy types, using the spec- 
tral energy distributions (SEDs) from the Galaxy Isochrone 
Synthesis Spectral Evolution Library model ( GISSEL96, see 
lBruzual&Charlotlll993t lLeifhereretaI1ll996l) . The code de- 
rives a matrix, N[m(F),z,t], containing the number of galax- 
ies at each redshift z (per square degree and apparent magni- 
tude interval) that contributes to the galaxy number counts at 
apparent magnitude m, for a given band F, and for the mor- 
phological type t. K- and e-corrections (evolutive corrections) 
corresp onding to each galaxy typ e t at each redshift z are con- 
sidered ( Yoshii & Takahara 1988), as well as dust-extinction ef- 
fects (iBruzual et al.ll988HWangll99lL see § 2.6 by G98 for more 
details). The original code implementation considers that galax- 
ies do not experience any morphological evolution, so it just es- 
timates the luminosity evolution backwards-in-time of a given 
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present-day morphological type, according to the SFH assumed 
for this type. 

G98 also inserted a procedure to account for number evo- 
lution in his model. Its merging procedure simply removes a 
given fraction of galaxies at each redshift from those galaxy 
types that the user considers to be affected by mergers. In this 
regard, the procedure is undoing the remnant galaxies into their 
original progenitors (we will refer to this as reversed merging). 
This procedure has some problems in the counting method (as 
shown in 92.2.11 ) and the merging procedures are too simplistic 
to account realistically for mergers. In particular, the assump- 
tions that mergers can only occur between galaxies of the same 
galaxy type, that the remnant must also be of this type, and that 
the merger is instantaneous are clearly contrary to observations 
(see assumption #2 in 92.3.21 ). Moreover, these procedures do 
not account for the cumulative effects of mergers on the LFs ei- 
ther. So, we have improved the original code to account more 
realistically for the effects of major mergers in galaxy evolution 
(see 92.31 ). For more details on the original code, the reader is 
referred to G98. 

2.2. General improvements to the model 

We have modified the original NCMOD code to derive the LF of 
each galaxy type at each redshift, <p[M(F), z, t\. In our improve- 
ments to the code, we have studied the morphological evolu- 
tion only driven by the major mergers strictly-reported by ob- 
servational studies. We are just interested in testing if the major 
mergers reported by observations can really provide an explana- 
tion to the observed evolution of LFs, or if additional processes 
are required to explain it. As a major merger basically converts 
the two progenitor merging galaxies into an E-SOa system (see 
assumption #2 in 92.3.2I ), our code traces the effects of major 
mergers on galaxy evolution backwards-in-time just decompos- 
ing present-day early-type galaxies (ETGs) into their progenitors 
(which are of different morphological type). So, the only mor- 
phological transformation between types allowed in our code is 
through reversed major mergers, i.e., from a ETG to its progen- 
itors. 

This means that we have not accounted for other processes 
that can lead to a transformation between types (such as secu- 
lar evolution through bars or minor mergers). If their contribu- 
tion to the buildup of massive ETGs were really significant, our 
model should be unable of reproducing the observed evolution 
of the galaxy LFs just accounting for the effects of major merg- 
ers. However, observations pose that these processes have been 
more relevant in the formation of intermediate -mass disks of 
types later than Sa than in the bui ldup of ETGs (lAzzollini et alJ 
2009; Lopez-Saniu an et al.l 12010). As the massive end of the 
galaxy LFs is controlled by ETGs up to z ~ 1, the assumption 
that Sa's-Irr's do not experience any morphological transforma- 
tion backwards-in-time does not affect to the results derived for 
massive systems. We comment below the improvements imple- 
mented to the original code in detail. 

2.2.1. Counting method 

Most of the authors using the backwards-in-time technique 
assume that the number density of galaxies of type t at z, 
with absolute magnitude in the F-band M(F), is equal to the 
number density of this galaxy type at z — exhibiting the 
same abs olute magnitude, i.e., <p\M(F),z, t] - < j>\M(F),z = 
0, t] (see Mas lankall 19841 lYoshii & Takaharall 19881: IWandll99U 



IKitzbichler & Whitell2006l and eq. 3 in G98). However, this as- 
sumption does not account for the effects of the L-evolution of 
galaxies. In fact, the SED of a certain galaxy type at redshift z 
has changed with respect to that at z = 0, because of its SFH. 
Therefore, galaxies exhibiting M(F) at z do not exhibit the same 
M(F) at z = 0, unless no L-evolution is assumed, and thus: 
<f>[M(F),z, t] * 4>[M{F),z = 0, t] in general. 

The relation between M(F) at redshift z and the absolute 
magnitude that a galaxy of type t exhibit at z = 0, M(F, z = 0), 
is given by the following expression: 

M(F)-M(F, Z = 0) = -2.51o g ff£|^, (2) 

where F(A) is the filter throughput, and fA[tc(z)',t] and 
//i[ f c(0); t] are the SEDs of a galaxy of type t at redshift z and at 
present, respectively. So, the galaxies observed at z with M(F) 
would be those that appear in the local LF with absolute magni- 
tude M(F, z — 0) at z — (according to the previous equation), 
and not those with M(F). This means that, in case of PLE: 

cf>[M(F),z,t] = 4>[M(F,z = 0),z = 0,f]. (3) 

If we also have number evolution, the previous equation is 
not valid, as galaxies from the original LFs at z = are dis- 
sapearing backwards-in-time through mergers. In this case, the 
cumulative effects of number evolution from z = to z must be 
considered: 

cp[M(F),z,t] = 4> m [M(F,z = 0),z = -» z,t], (4) 

where <p m [M{F , z — 0), z — — » z, t] represents the number den- 
sity of galaxies of type t still existing at z, having into account the 
backwards-in-time effects of major mergers since z — up to z 
in this population. Tests performed to estimate the effects of the 
different counting methods have shown that the old procedure 
underestimates the galaxy number counts at faint magnitudes as 
compared to our method, by a factor of ~ 25 at m(B) = 25 mag, 
and twice more at m(K) = 16 - 18 mag. Effects start to be sig- 
nificant at z > 0.8 in all bands (from U to K). 

2.2.2. Local luminosity functions by morphological type 

Galaxies within t he same morphological class also exhibit simi- 
lar colors. In fact. lKauffmann et al.l (120041) show that galaxy mor- 
phology is strongly related to star formation (SF). Therefore, we 
have decided to follow the evolution of the different morpholog- 
ical types existing at z ~ backwards-in-time, assuming a priori 
that all the galaxies within each type share assembly processes 
and SFHs. 

The code requires as input the local LFs of each galaxy 
type. Two studies on the local LFs given by morphological types 
were available: 1) the LFs derived by iNakamura et ail |2003, 
N03 from now) in the r*-band, provided in four coarse morpho- 
logical classes (E-SO, SOa-Sb , Sbc-Sd, Im); and 2) the LFs by 
ICuesta-Bolao & Sernal (120031 C03 hereandafter) in the B z filter, 
for eleven thin morphological classes (E, SO, SOa, Sa, Sab, Sb, 
Sbc, Sc, Scd, Sd, Im). None of these two parametrizations fit- 
ted our requirements, as the last one gives better predictions on 
those observables associated to star-forming systems, whereas 
the former was better tracer of the mass evolution. Moreover, 
the SOa-Sb class by N03 is to o coarse as to be represent ed with 
the same local SED and SFH (lEliche-Moral et al.(l2006bh . 

So, we have decomposed the r*-band LFs by N03 into thin- 
ner classes, using the information provided by the B-band LFs by 
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M(r*) (mag) 

Fig. 1. Decomposition of the local LFs by N03 in r* into thinner 
morphological types, using the local LFs in B by C03. The LFs 
plotted correspond to the following morphological classes: E- 
SOa, Sa-Sab, Sb-Sbc, Sc-Scd, and Sd-Irr (consult the legend in 
the figure). Solid line: Total LF resulting from the performed 
decomposition. Thick, dashed line: Original total LF by N03 in 
r*-band. [A color version of this plot is available at the electronic 
edition] . 



C03, and trying to group more accordingly galaxies with simi- 
lar spectrophotometric and morphological properties (E-SOa, Sa- 
Sab, Sb-Sbc, Sc-Scd, and Sd-Irr). The method is explained in 
detail in AppendixlAl FigureQ] shows the resulting local LFs in 
r*, as derived using this procedure. Notice that the LF obtained 
for the E-SOa class exhibits an increasing slope at the faint-end. 
This does not contradict studies repo rting a decreasing fa int-end 
slope (a > -1.0) for the E-SO class (de Lapparent 2003), as our 
class also includes the SOa's. Schechter fits have been performed 
to these decomposed LFs for each morphological type, obtain- 
ing the parameters and^- 2 residuals listed in TableQ] These LFs 
in the r*-band have been considered as the input local LFs in the 
model. They trace the SF better than the original LFs by N03, 
with the advantage of a better correspondence with the stellar 
mass, as compared to the original LFs in the B-band. Several 
tests to check the robustness of this decomposition are com- 
mented in AppendixlAl 

2.2.3. Star formation histories 

The SFH of each galaxy type in the model is fixed through the 
star formation rate (SFR) assigned to each type and the redshift 
at which the SF activity in each galaxy class starts (zt>). We 
have used an exponential fun ction to represent the SFRs of spi- 
ral galaxies dKennicutd ll989). a single stellar population (SSP) 
for ETGs, and a constant SFR for Sd-Irr's. Different e-folding 
timescales have been considered for spirals, depending on their 
types and according to observations: t , = 4Gyr for Sa -Sab, 
7 Gyr for Sb-Sbc, and 9 Gyr for Sc-Scd (iGallazzi et alJl2005h . 



of mETGs at z < 1 with mass-downsizing 



Table 1. Schechter Parameters of the decomposed LFs a in r* 



Type 


f (0.01Mpc- J ) 


M'(r') 


a 




(1) 


(2) 


(3) 


(4) 


(5) 


E-SOa 


0.20±0.22 


-21.51±0.20 


-1.0 b 


7.8e-ll 


Sa-Sab 


0.3±0.4 


-20.62±0.03 


-0.49±0.14 


1.3e-10 


Sb-Sbc 


0.169±0.007 


-21.27±0.04 


-1.13+0.03 


8.7e-10 


Sc-Scd 


0.1±0.4 


-21.06±0.04 


-0.69±0.15 


1.4e-10 


Sd-Irr 


0.037±0.020 


-20.94±0.07 


-1.90±0.03 


2.0e-9 



a Provided for h = Z/o/100 = 1. Magnitudes in the Vega system. 
b Fixed a value (see AppendixfA). 



Table 2. SFRs, metallicities, and internal dust extinction per 
galaxy type 



Galaxy Type 


SFR 




Z b 


T imU t*(B, @R 25 /2f 


(1) 


(2) 


(3) 


(4) 


(5) 


E-SOa 


SSP 




0.02 


0.10 


Sa-Sab 


exponential 


4 


0.02 


0.55 


Sb-Sbc 


exponential 


7 


0.008 


0.70 


Sc-Scd 


exponential 


9 


0.004 


0.37 


Sd-Irr 


constant 




0.004 


0.09 



a e-folding timescale in Gyr. 
b Metallicity. 

c Optical depth in B due to internal dust extinction (M09). See the 
text for more details. 



We have assumed metallicities that decrease from early 
to late types, in ag reement with the observational results by 
Gall azzi et al.l(l2005l) . The final parametrizations of the SFR con- 
sidered for each galaxy type in the model are listed in Table|2] 
In fact, these values are qui t e standard in population synt hesis 
models dGallazzi et al.ll2005t lTosill2009l: iMaier et alj l2009). and 
similar to those used in lEliche-Moral et all d2006bf) . The ages 
and metallicities considered in the model are coherent with the 
median values of thes e magnitudes derived for SDSS galaxies in 
different mass ranges (Gallazzi et al. 2005h. 

For all galaxy types, we assume that the SF starts at 
zs.t _ 3, where the SFH of the Universe exhibits its peak 
dHopkins & Beacom 2006). Notice that this redshift represents 
the formation redshift of the bulk of the SPs for all the morpho- 
logical types. This is different to the redshift at which the galaxy 
is assembled, as stars can be formed in different systems to those 
where they end up at z = 0, because of merging. 

To simulate the spectral evolution of the composite SPs we 
have updated the library of evolutionary SP synthesis models of 
the original code to GALA XEV, the Isochrone synthesis code by 
iBruzual & Charlotl d2003l B&C03 hereandafter). Padova 1994 
evolutionary trac ks and the Chabrier initial mass function have 
been considered dChabrieril200l I2002L l2003n . 

2.2.4. Dust extinction 

G98 considered all the galaxy types to be equally extinguished, 
using a B-band optical depth T"d U st ; star(B) = 0.2 for local L* galax- 
ies of all types at 4500 A. This assumption is a sim plistic ap- 
proximation to reality dBuat et al.ll2002t iBoselli et alj|2002l) . So, 
we have inserted of a different Td us t,star(B) for each morpholog- 
ical type. We have assumed the average values of dust extinc- 
tio n derived from the observa tional extinction profiles obtained 
by|M 

unoz-Mateos et al. ( 2009, M09 hereafter) for each morpho- 
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logical type. We have adopted the extinction value at its isopho- 
tal radius /?2s/2 as the representative total internal extinction of 
that galaxy type (see Table|2|. The assumed Td us t,star(B) values 
are consistent with those derived by other authors within the typ- 
ical o bservational errors (+0.1, see iBuat & Xulfl996llTrewhenal 
1998). 



2.3. Merging 

In the model, we are going to consider the effects of major 
mergers on the galaxy populations, but not those of minor 
mergers. Although th eir role on the galaxy evolution is poorly 
under s tood at present dFerreiro & Pastoriza 2004; Gimeno ~et al.l 



2007; Martmez-Delgado et al. 2009), we know that major 
mergers drive much more dramatic structural changes, vio- 
lent starbursts, a nd higher mass incr e ments in a galaxy than 
minor mergers dBarnes & Hernquisa 1991 ; [B endo & Barnes 
2000; Naab & Burked 12003b iGonzalez-Garcfa & Balcellsl 
20051: lEliche-Moral et all l2006at IGonzalez-Garcfa et al.l 12006 



Bournaud et al. I2007I) . Even considering that minor mergers 



have probably been more numero us than major mergers (at 
least, by a factor of ~ 2, s e e e.g.. lEliche-Moraletal.ll2 006a; 



Hetznecker & Burkert 2006; Bournaud et al. 2007; Cox et al. 
2008bt iKeell 120081: IStewart et alJ 120081: iJogee et al.l 120091: 
Kavirai et alJl2009bi i Naab et all 120091: IStewart et alJ l2009bTa[ 
Lopez-Saniuan et all 120101) . observations indicate that their 
contribution seems to have been signifi cant only in low- 



mass systems (wi t h MJAU < 10 10 , see iBundv et al.l [2009; 



lOesch et alJ l2009t lLopez-Saniuan et alJ |2010|). Therefore, 
can neglect their effects in the model. 

2.3.1 . Evolutionary stages in a major merger 

Depending on the gas content of the progenitors in a major 
merger, their colors can redden noticeably due to strong dust 
extinc tion during in termediate phases of th e merger ( Darg et al.l 
I2010I) . According to lLotz et all d2008bl) and lCox et al l d2008al) . a 
wet major merger reddens noticeably due to strong dust extinc- 
tion during two epochs of strong merger-induced starbursts: 1) 
in the pre-merger stage or tidal phase (enclosed between the first 
pericenter passage and the moment at which both progenitors 
have merged into one body), and 2) in the merging-nuclei stage 
or nuclear coalescence phase (during which both nuclei merge 
generating a unique remnant). In both phases, the merging galax- 
ies also exhibit high morphological distorti ons (iLeithereil 2000. 
200 U iBell et al.l l2006bi IConselicd l2006bt lEmonts et aljbooa 
Hernandez-Toledo et alJ|2006l; ICox et al.ll2q08at Idi Matteo et all 
20081: lLotz et al I l2008albl: IConselicd l2009t iGeach etalJ 12009: 
Kim et all 120091) . During the last stage of the merger (post- 



merger or late remnant phase), the SF in the remnant is com- 
pletely quenched by the galactic superwinds associated to the 
recently-formed stars or by AGN feedback, the SPs in the rem- 
nant fade passively, and the system experiences a strong relax- 
ation (see iHeckman et apll999t lEmonts et al. 2006: ICox et alJ 
12004 12008at lLotz et alj|2008bl: ISchawinski et all 120071. 120101) . 
If gas ejection from the remnant is efficient enough and gas in- 
fall is negligible (as at it seems to be usual at z < 1), the remnant 
evolves pas sively into an ETG in T3 ~ 0.5-1 Gyr after the SF 
quenching dNaab & Burkert|l2003l: ISchweizei|l2005t INaab et al . " 
2006; Robertson et afll2006: Ems ellem et alJl2007t iRossa et al 



20071: ISchweizer & Seitzerll2007l) . 

A mixed major merger presents the same structural stages 
as a wet merger, although the gas content and SF in- 



volved are lower. As it involves a gas-rich, spiral pro- 
genitor, appreciable morphological distortions make it de- 
tectable to asym metry ind i ces stu dies during the nuclei- 
merging phase dLotz et alJ l2008bl) . The starbursts driven 
by the int eraction can al so re dden noticeably the rem- 
nant center dHearn & Lambl 120011; IStatler & Mc Namaral 120021: 
iHernandez-Toledo et al.ll2006HJohansson et al.ll2009l) . 

A dry major merger imprints negligible morphological 
distortions and entails very little SF, and thus, the remnant 
and the progenitors can be considered typical ETGs during 
the w hole process (jNaab et alJ 120061: Hern andez-Toled o et al.l 
2006 ; IBell et al ll2006ab iRines et alj|2007t iDonovan et al.l l2007: 
Feld mann et alJl2008l) . 

The nominal timescales assumed for the pre-merger and 
merging-nuclei phases correspond to the ave rage values derived 
from the whole range of values reported by lLotz et al.l (|2008b) 
and ICox et"all d2008al) : t\ ~ 0.7 and r 2 ~ 1 Gyr, respectively. 
However, we have also accounted for the uncertainties inserted 
in the model due to these timescales, considering the extreme 
values reported by these authors for each timescale (see $3] f° r 
more details). 



2.3.2. Assumptions on the merging events 

Concerning to major mergers, we have assumed in the model 
that: 

1. A major merger is a collision between two galaxies, 
involving stellar mass ratios ranging between 1:1 to 1:3. — 
This is supported by the observational fact that major merg- 
ers involving more than two g alaxies simultaneousl y are ex- 
tremely rare in the Universe dVaisanen et al.l 120081) . and by 
recent studies stating that merger fractions determined us- 
ing concentration-asymmetry-clumpiness (CAS) methods are 
only sensitive to major mergers involvin g stellar mass ra- 
tios between 1:1 and 1: 3 dConselicdl2006allbl : lLotz e7al1 f2004. 
2008b). In hierarchical models, une qual-mass mergers are more 
frequent than e qual-mass mergers dKhochfar & Burkertl l200ll 
Khochfar & Silk 2006), but as the distribution of mass ratios in 
real mergers for a given mass is unknown, we have assigned a 
mass ratio 0.3 < mr < 1 to each merger event according to a 
random distribution. 

In order to establish a direct relation between stellar masses 
and luminosities in galaxies, we have performed all the computa- 
tions in a NIR band (K). This assures that the stellar mass At* of 
a galaxy is nearly proportional to the luminosity L in that band, 
whatever the morphological type of the galaxy is. Neglecting the 
effects of merger-induced SF in first approximation, the absolute 
/f-band magnitudes of both progenitors can be derived as: 



M X (K)-M(K) = -2.5 log[l /(mr+ 1)] 
M 2 (K) - M(K) = -2.5 log[(mr/(mr + 1)], 



(5) 



being M^K), M 2 (K), and M(K) the absolute K-band mag- 
nitudes of the two progenitors and of the resulting ETG, 
respectively. 

2. A major merger gives place to an ETG, independently 
of the morphological types of the progenitors. — There are nu- 
merous observational results that point to a major merger origin 
for most elliptical gal a xies dSchweizerlj2005:jRossa et aL 2007 



Feld mann et ail 1200 8; Poggi anti et alj|2009at iFalkenber g et al.l 
2009). Although numerical simulations show that the final rem- 
nant structure depends strongly on the mass ratio, gas content, 
and orbit of the encounter, they state that a major merger be- 
tween galaxies with typical gas contents basically gives place 
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Fig. 2. Relative fractions of wet, mixed, and dry mergers as a 
function of redshift reported by L08 (squares, triangles, and di- 
amonds, respectively). Lines: Linear extrapolations of the values 
at z — and z = 1. Consult the legend in the figure. Vertical 
dashed lines: Transition redshifts at which the fraction of one 
of the merger types becomes null, according to the extrapola- 
tions. [A color version of this plot is available at the electronic 
edition] . 



to an E-SOa dBarnes & Hernquisll [1996] iBournaud et al] 120041 



120051: iNaab & Burkeril2003t iNaab et all2006l) 

This assumption does not contradict t he di sk rebuild- 
ing scenario proposed by [Hammer et al.l (120051) (see also 
ISpringel & Hernquistl 120051: Hopkins et ail 120091) . According 
to it, a large disk could be rebuilt after a major merger in 
~ 2-3 Gyr, in the case that the progenitor galaxies contain a 
large gas reservoir (amounting to ~ 50% of their masses, see 
Hammer et al. 2009a: iPuech et al.l l2009a). This s cenario is sup- 
ported by observational examples at 0.6 < z < 1 (iHammer et aF 



2005: iGovernato et al.ll2009t IHammer et alll2009bt IPuech et al. 
2009b. lYang' et al.l 120091) and by numerical and cosmological 
simu lations (see nume rous references in the introduction 
by IStewart et al.l 12009 al) . However, a close inspection to the 
sufficiently-relaxed observational examples shows that those 
cases related with major mergers are basically SO-SOa galaxies, 
whereas those exhibiting lower bulge-to-disk ratios (and thus, 
not being ETGs accor ding to our definition) are associated to 
minor mergers inste ad (Elich e-Moral et ai1l2006al: ICassata et al.1 
2008; iPoggianti etaP l2009al) . Therefore, the disk rebuilding 
does not conflict with our main hy pothesis at z, < 1 (see an 
interesting discussion on the topic in lOesch et al.ll2009h . 

3. Major merger fractions vary as a function of redshift. — 
Observational errors usually associated to typical CAS methods 
lead to a p ossible overestima t ion of merger frac tions by a factor 
up to ~ 2 (lLotz et alj|2008at iJogee et alll2009l). In order to ac- 
count for these uncertainties, Lopez-Saniuan et al. (2008) have 
developed a maximum likelihood (ML) method that considers 
the biasses linked to the redshift degradation of the morphologi- 
cal information. This tech nique provides lower merger fractions 
than typical CAS methods dLopez-Saniuan et alJ2 009a b. LSJ09 
hereandafter), quite in agreemen t with visual and pair count- 
ing estimates (Jogeeetal. 2009). Therefore, we have adopted 
LSJ09 merger fractions up to z ~ 1 for the model, keeping in 
mind that they are limited to galaxies brighter that M(B) - -20 
(Vega magnitudes). This corresponds roughly to galaxies with 



M < M* for z < 1 . According to these authors, the merger frac- 
tion evolves as a power-law at z < 1 as: 



(6) 



where m = 1.8 + 0.5 and / m (0) = 0.014 + 0.003. The un- 
certainties of these parameters have been considered for error 
assessment (see $3). The extrapolation of LSJ09 trends at z > 1 
provid es similar results to those derived by Con selice et al.l 
(2008) at 1 < z < 1.5. However, the high observational errors 
of merger fractions at z > 1 make the model results at these 
redshifts uncertain. 

4. Merger fractions in the B-band trace reliably the number 
of merger events at z 55 1 ■ — Some studies report that merger 
fractions derived in NIR are lower than those estimated from 
optical bands, probably because the star-forming features of 
galaxies seem more asymmetrical as we move to shorter rest- 
frame w avelengths dBundv et al.ll2.004l : iHuertas-Companv et al 
2009b; Rawat et al.l 120091) . However, lLopez-Saniuan et al 



(2009a) find that the morphological merger fractions at z = 0.6 
in the K and B bands are compatible within the errors. 

5. Merger fractions as a function of redshift for different 
merger types. — Mergers involving gas-rich (late-type) galaxies 
dominate at z > 0.8, whereas mergers involving gas-p oor (early- 
type) galaxies are more f r equent at recent epochs ( Tran et al. 



type) galaxies are more trequent at recent epocns ( Iran et al. 
2005: Ivan Dokkuml 120051 iLin et al] 12008b lLotz et al] 12008a: 



de Ravel et all l2009UBundv et al.l 12009). Percentages of wet 
mixed, and dry mergers have been estimated up to z ~ 1 by 
ILin et al] (120081 hereandafter L08). We will consider their 
relative fractions in order to estimate the net number of wet, 
mixed, and dry mergers taking place at each redshift. Figure 
[2] shows the linear extrapolations in redshift derived from L08 
values at z ~ 0.1 and z ~ 1 that have been considered in the 
model. 

6. Merger fractions by LSJ09 are only sensitive to the 
merging-nuclei stage of gas-rich mergers. — Observations have 
stated that CAS methods are not efficient detecting dry mergers 
van Dokkuml 120051: [Bell et al] l2006bt iHernandez-Toledo et all 



120061: iDonovan et al.l 12007 



Feldman n et al.l |2008). Therefore, 
we will assume that the merger fractions by LSJ09 only account 
for wet and mixed processes. Under this assumption, eq.Q] be- 
comes: 



0m,w(z) + 0m,M(z) = /mfe) ■ 0total(z), 



(7) 



where f m represents the merger fraction by LSJ09 and m ,w and 
0m,M the number density of wet and mixed mergers detected at z. 
The relative fractions of wet, mixed, and dry mergers at z are es- 
timated as the number of mergers of each type detected at that z, 
divided by the total number of mergers at that z (L08). Therefore, 
the relative fractions of dry, mixed, and wet major mergers can 
also be computed as: 

/d(z) = 0m,D(z)/0m(z), 

/m(z) = 0m,M(z)/0m(z), (8) 
/w(z) = 0m,w(z)/0m(z), 

where (f> m (z) and m ,D(z) represent the number density of all and 
dry major mergers at z, respectively. The total number of merger 
events detected by L08 also contain dry mergers, and thus: 

(pm(z) = 0m,D(z) + 0m,W+M(z), (9) 

with m ,w+M accounting for the sum of m> w and m> M- 
Therefore, dividing eq.|7]by (f> m (z) and replacing the definitions 
given in eqs.[8]|9] we can estimate the total number density of 
major mergers at redshift z through: 
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a(z) = /m(z) ' <f>total(z)/[fu(z) + /W(z)L 



(10) 



as well as the number density of dry, mixed, and wet mergers at 
z through eqs.[8] 

CAS methods detect me rgers in the immediate pre-merger 
and merging-nuclei phases dLe Fe vre et al. 2000; Patton et al. 
2002tlConselice et aljpOOllBell et alJl2006bHConselicel l2006b: 
De Propris et alj|2007t Ide Ravel et all 120091) ! Nevertheless, the 
visual inspection of the merger samples by LSJ09 indicates 
that the mergers correspond basically to merging-nuclei stages 
of gas-rich encounters (probably due to the image degradation 
performed by these authors). So, we have assumed that all the 
mergers included in LSJ09 merger fractions have been detected 
at this evolutionary stage. 

7. Merger fractions depend on the stellar mass. — The 
merger rate increases as the galaxy mass increases at 
Z < 0.6, whereas at higher redshifts the merger rate does 
not depend on mass for galaxies with A(* > 10 n A1o 



(Carlberg et al. 


1999, 2000; Parkinson et al. 


2008; Genel et al. 


2009 


; Bundv et al.1 2009; Conselice et alj 20091 de Ravel et al. 


2009 


; Lopez-Saniuan et al. 2009a; Domingue et alJ 2009). 



As the model is limited to bright galaxies (L > L*), we 
will distribute the computed number of merger events at a 
given redshift according to these trends. This means that the 
fractional increment in stellar mass at each mass bin is the 
same for all masses in the model at z > 0.6 by assumption 
(AAL/At* = constant, VA1*), i n agreement with observa- 
tions and semianal ytical models dDe Luc ia & Blaizot 2007; 
iFontanot et al]|2009h . 

8. Starbursts in gas-rich major mergers. — Observations and 
simulations have posed that mergers between gas-rich galaxies 
induce violent starburst s in the remnant, widely-spread along 
the whole galaxy dJoseph & Wrighl 1985b iMihos & Hernquisl 
19961: IZezas et alJ 2003t iKannappan et all 120041: iBridge et alJ 
20071 Ijogee et alJ 120091: iKim et alJ 120091: iDarg et alj l2010). 
The efficiency of mergers in triggering strong starbursts is high 
enough as to consume ~ 50-75% of the original gas content 
dCox et alJl2004 120061) . Therefore, we have assumed that: 1) 
each major merger involving a gas-rich progenitor must have 
converted a certain percentage of its original gas content into 
stars, and 2) the starburst affects the color of the whole galaxy. 
Assuming that the typi cal gas content in disks is ~ 10% of 
their baryonic masses (Robertsl [19751 Young & Scoviilelll982l: 
IScoville & Yound ll983). we have tested how the efficiency of 
the merger-induced starbursts can affect to the predicted LFs 
(see 92.51 l. We have considered that the increment in stellar 
mass induced by the merger is concentrated at the pre-merger 
phase. Assuming that a fraction / (0 < / < 0.1) of the original 
gas content of each disk progenitor has transformed into stellar 
mass in the remnant ETG, a stellar mass Al* (1 — f) (instead 
of a mass Al*) must be distributed into the original gas-rich 
progenitors, affecting to eq.[5] 

9. DSF phases in gas-rich major mergers. — The strong 
starbursts induced by major mergers are frequently accompanied 
by large amounts of dust, that extinct and redden noticeably 
the galaxy spect r um dKeel & Wulll995t iBell & de Jondfe OQ 1 ; 
Kauffman netai] 120031: iRoche et alJ 120061: iPoggianti et all 
2009ah . In the local Universe, gas-rich major mergers ex- 
hibit higher dust contents than normal star-forming disks 
dBekki & Shioval 120001: IPollack et al J 120071: Ilmanishill2009l) . In 
fact, colors of local ULIRGs (which are basic ally gas-rich major 
merg e rs) are as red as those of local ETGs (Surace & Sanders 
2000; Sur ace et alfe OQO: Arr ibas et all 120041: lllbert et al.l l2006: 



IDarg et alJ 1201 Oh . At high redshift, major mergers exhibit 
even higher amou nts of SF and dust extinction than their 
local counterparts dCalzettH 1 1 999b lAlonso-Herrero et ail [2004; 
ICortese et all 120061: llono et alJ l2009h . Consequently. ETGs 
and DSFs at high redshift exhibit similar SEDs, except for 
some abs orption and emission lines present in star-forming 
systems dLiu et alJl2000b iPierre et alJl200llls"trateva et al lfeOQ 1 ; 
I Cimatti et al.l|2002l;IWeiner et alll2005l;lvan Dokkum et al| 2006; 
Ischurer et al] l2009t) . In fact. ISchawinski et all d2007M 20 1 0) 
have proved observationally the existence of an evolutionary 
sequence in local ETGs from a merger-driven starburst phase 
to quiescence, that lasts ~ 1 Gyr. The transition between both 
phases is driven by an AGN phase that peaks at ~ 0.5 Gyr after 
the starting of the starburst phase, responsible of the suppression 
of the SF in these systems. These authors report colors of ETGs 
in the AGN phase quite similar to those exhibited by quiescent 
ETGs on the red sequence. Therefore, we will assume that 
gas-rich mergers undergoing transitory DSF phases (during the 
pre-merger and merging-nuclei phases) exhibit the same SED 
as ETGs at the same redshift. This assumption is reasonable 
in global terms at least for rest-frame wavelengths longer than 
the 4000 A break feature. Moreover, we will show that DSFs 
become numerically-relevant at z > 0.7 in the model (see §0, 
epochs at which real DSFs exhibit the reddest colors, making 
our assumption more accurate. 

10. The gas-rich progenitors in a wet or mixed major merger 
are usually spirals of very late Hubble type. — This assump- 
tion is supported by observations posing that gas-rich galax- 
ies involved in wet or mixed mergers tend to be of Sc-Irr type 
(De Propris et al. 2005), and by the fact that the relative fraction 
of early-type galaxies in high-density environments seems to 
have increased since z ~ 1 at the expense of that of ver y late-type 
galax ies (Sc-Irr's), instead of intermediate spirals dTasca et all 
2009). So, we assume that, if an ETG was built up through a 
mixed merger, its gas rich progenitor was a Sc-Irr galaxy (the an- 
other one was a gas-poor galaxy), and if the ETG was assembled 
through a wet merger, both gas-rich progenitors were Sc-Irr's. 

2.4. Computing the LF evolution 



In our model, the backwards-in-time evolution of the local LFs 
is computed assuming that only ETGs are susceptible of com- 
ing from a major merger. In order to derive the evolved LFs at a 
given redshift z, we carry out a procedure analogous to the orig- 
inal G98 method ( 92.1I >. but using the correct counting method 
( 92.2. U and all the improvements commented in 92.21 

We assume that the number density of E-SOa's that are 
formed at each redshift z through major mergers is equal to the 
number density of major mergers that are observed at z (assump- 
tion #2 in 92.3.2b . The number density of the assembled ETGs at 
each redshift is estimated using LSJ09 merger fractions (eqs.|6] 
and [TOt. The corresponding number densities of ETGs formed 
through wet, mixed, and dry major mergers at each redshift are 
estimated through eq[8] The number of ETGs that must be re- 
moved from their LF at each luminosity bin (and at each redshift) 
is estimated assuming the distribution of mergers with luminos- 
ity stated in assumption #7 of 92.3.21 As the merger fractions by 
LSJ09 are detecting remnants before their post-merger phases, 
we can override this phase in the model (assumptions #3-6 in 
9232). 

We have decomposed each ETG of luminosity L coming 
from a major merger into two progenitors, with luminosities ac- 
cording to equations derived in 92. 3 .21 (assumptions #1 and #8). 
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In reversed merging, these recently-assembled ETGs must dis- 
appear from the cosmic scenario backwards-in-time, but firstly 
they come into the nuclei-merging phase during r 2 Gyrs, split 
into their progenitors in the pre-merger phase for t\ Gyrs, and 
then become two typical progenitor galaxies (ETGs or blue 
disks) after this phase (as described in 32.3. 11 1. This procedure 
has into account that the conversion of two merging progeni- 
tors into an ETG is not instantaneous, and the colors of galaxies 
at each phase depending on the merger type. The product of a 
reversed merger in our back-in-time technique results in the de- 
composition of the mETG assembled through this merger into 
their two progenitor galaxies (ETGs or gas-rich disks, depend- 
ing on the merger type), after a period of t 2 + T\ Gyr. Gas-rich 
progenitors have been considered to be Sc-Irr galaxies (assump- 
tion #10 in 32.3.21 ). Galaxies undergoing a DSF phase during a 
gas-rich major merger (in the pre-merger and in the merging- 
nuclei phases) are considered to exhibit the same SED as ETGs 
at each redshift (non-corrected from dust extinction, see assump- 
tion #9 in 32.3.2b . Notice that, as the model traces the evolu- 
tion backwards-in-time, the fate of the most massive galaxies is 
not affected by the evolution of the less massive systems, as all 
galaxies are continuously losing stellar mass backwards in time. 

The timescales adopted for the DSF phases of wet and mixed 
major mergers control the "life period" of the transitory DSFs 
before becoming typical late-type disks, and hence, affect to 
the LF of these systems. However, it is worth noting that these 
timescales (ri ~ 0.7 Gyr for the pre-merger stage and t 2 ~ 1 Gyr 
for the merging-nuclei phase) do not affect the net assembly 
rate of mETGs at each redshift, but they just delay more or 
less the transformation of the progenitor galaxies in transitory 
merging stages into "typical" systems. The number of mETGs 
dissapearing at each redshift depends only on the merger frac- 
tions, the relative percentages of each merger type, and the total 
galaxy LF at each redshift, but not on any timescale (see eq.fTOb. 
As observational merger fractions are independent on author's 
hypothese s on the timescale of a merger (contrary to merger 
rates, see iLe Fevre et alj 120001 : iBell et all l2006bt iJogee et ail 
2009, LSJ09), the number of mETGs assembled through major 
mergers estimated by our procedure is independent on the se- 
lected timescales of the pre-merger and merging-nuclei phases 
(at which the built up ETG does not exist, but their merging pro- 
genitors). 

2.5. Model uncertainties and error assessment 

The values indicated in § 32.1112.31 for parameters such as the 
SFH per morphological type, the merger fractions, or the merg- 
ing timescales are average values derived from observations and 
simulations. We have tested how errors in the nominal model 
parameters affect the model results. 

i- SFHs for each morphological type ( 32.2.31 ). — 
Eliche-M oral et al.l d2006bi) showed that, although galaxy 
number counts predicted with NCMOD are not sensitive to the 
SFH selected for each galaxy type, the predicted color distribu- 
tions depend noticeably on them. We have tested how the LFs 
predicted by our model change when the standard parametriza- 
tions assumed for the SFH of each galaxy type are modified. 
Our tests confirm that the resulting LFs are quite robust against 
small changes of these parameters, provided that the input 
values are nearly realistic. Moreover, the morphological types 
used in the present mode l group galaxies w i th mor e similar 
spectral properties than in lEliche-Moral et all d2006bl) . making 



more realistic the representation of all galaxies within a given 
type with the same SFH (see 32.2.2l and AppendixlAl. 

The SFH assigned to the Sa-Sab's gives place to slightly 
bluer colors for these galaxies than those exhibited by this 
population (see 33.21 ). If we consider a shorter timescale for 
its SFR, in order to make them redder, the fitting in NIR 
bands of their LFs (and even of galaxy number counts, see the 
forthcoming paper III of this series) improves, but the LFs of 
the total galactic population are a bit overestimated. This galaxy 
population has evolved probabl y through secular ev o lutionary 
process from late-type disks (lAthanassoulal 120081 : ICombes 
2008; Kormendy & Fisher 2008), processes that are not being 
considered in the present model. So we have fixed their SFH to 
the typical values reported by observations for this galaxy type. 

ii- Dust extinction for each type ( 32.2.41 ). — An accurate 
description of extinction is critical for comparing model predic- 
tions with data for bands from U through K. We found G98's 
assumption on dust extinction (Td ust (B) = 0.2 for all galaxy 
types) too conservative attending to studies of local disk galax- 
ies (see 32.2.41 ). We have tested that, for typical observational 
uncertainties of the optical depth (ATd us t,B ~ ±0.2 in the values 
of Tabled, results do not change appreciably. Higher dust 
extinctions in late-type galaxies lead to an underestimation of 
the LF of blue galaxies at all redshifts up to z ~ 1 . However, 
results in ETGs are less sensitive to dust extinction assumptions. 
In fact, they could tolerate dust extinction values beyond the 
above ATd us t,B ~ ±0.2, although the LF at 0.8 < z < 1 is 
worst-reproduced in this case. 

Hi- Assumed colors for DSFs ( 32.3.21 ). — We have assumed 
the same SED for representing DSFs and ETGs at each 
redshift, in order to simplify the model (see assumption #9 
in 32.3.21 ). This assumption implies that we are underesti- 
mating the UV emission of the DSF spectrum. So, we have 
not considered bands bluer than rest-frame B-band in our results. 

iv- Merger fractions ( 32.3.21 ). — These uncertainties have 
turned out to be the most relevant source of error in the 
model, so we show them together with results in 5j3] We have 
also implemented different the parametrizations of merger 
fractions apart from the one by LSJ0 9. In parti c ular, w e have 
tried the parametrization derived by IConselicel (1200 6a). This 
function provides a slightly better fit to LFs at 1 < z < 2 (by 
~ 10-20%). At z ~ 0.8, both parametrizations imply similar 
merger fractions. However, the scenario proposed by our model 
overestimates largely the number evolution of ETGs at z > 1, 
considering the major mergers reported by observations at z > 1 . 
In fact, we would require less than 20% of the major mergers 
derived from the observational merger fractions at z > 1 to 
reproduce the galaxy number counts by redshift bins and the 
LFs of red and blue galaxies at 1 < z < 1.5. Nevertheless, 
it seems that the properties of major mergers start to differ 
appreciably from their counterparts at lower redshifts at z > 1, 
and hence, assumptions adopted in the model are not fulfilled 
by real galaxies (see comments in 34. 31 ). Results are extremely 
sensitive to the adopted merger fractions (see §0 and the model 
uncertainty regions plotted in Figs.[3][T]). 

v- Merger timescales ( 32.3.71 ). — In order to account for 
the uncertainties in these parameters, we have considered 
the whole range of values reported by N -body simulation s 
(ti = 0.05- 1.7 Gyr and r 2 = 0.003-2 Gyr, see lCox et al.l l2008a). 
Although non-physical, we have also considered the possibility 
of "instantaneous" mergers (rj = t 2 = 0), in order to compare 
our results with those that would be obtained in this case, as 
assumed by several semi-analytical codes (see references in 
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92.2.11 ). We comment these uncertainties together with results in 
$3](see also the model uncertainty regions plotted in Figs.[3][T]l- 
vi- Merger-induced SF ( 92.3.21 ). — We have performed sev- 
eral tests changing the fraction of ETG stellar mass formed dur- 
ing the major merger that built it up. The typical increments of 
stellar mass stated by observations (betw een 0-10% of the to- 
tal stellar mass, see lCox et alJl2004T,l2006l) affects negligibly the 
final results. However, we have also tested how higher merger- 
induced star formation can affect to the results, finding that in- 
crements of stellar mass above ~ 40% of the total mass of the 
galaxies could affect the shape of the LFs of late-type galaxies 
at high redshift (see also item iii of 94.3b . 

3. Results 

In this section, we compare the evolution of the LFs predicted by 
the model for each galaxy population with observational data in 
different bands and according to different selection criteria (on 
color, on morphology). 

3.1. LF evolution of red and blue galaxies atz<l 

As pointed out by F07, studies focusing on red galaxies are im- 
plicitly assuming that rest-frame color is a good way of finding 
E-SO types at all redshifts. These authors claim that the iden- 
tification red galaxy-ETG "works well at low redshifts, where 
the contamination of the red sequence by DSFs wi th H ubble 
types later than SO is neg ligible (~ 10-15%, see iBell et al] 
2004a]; IWeiner et al.l l2005l) . However, contamination by non- 
spheroidal DSFs is larger at higher redshifts , amounting to 30% 
at z ~ 0-75, and may increase beyond t hat (Cimatti et al. 2002; 
Yan& Thompson 12003b IGilbank et al.1 l20o¥ " Moustakas et al. 
2004 ". [n fact, most of the authors selecting red galaxies admit 
that the similarity of the rest-frame colors of ETGs and DSFs 
makes di fficult to disentangle both galaxy populations at high 
redshifts (Strateva et al. 200 li lllbertet alJl200a ICirasuolo et al. 
120071; Franzett i et alj|2007l) . As the level of contamination by 
DSFs of a red galaxy sample is going to depend strongly on the 
color criteria used by each author, we are going to consider the 
two extreme situations in our model in order to compare with ob- 
servational data: 1) the case that the red galaxy sample includes 
exclusively ETGs, and 2) the case that all DSFs are also included 
into the red sample. 

3.1 .1 . fi-band LFs for red and blue galaxies up to z ~ 1 

In Figs.[3][5] we have compared the evolution of the rest-frame 
LFs in the B-band as predicted by the model (using our nominal 
parameters) with real data on red, blue, and total galaxy popu- 
lations. Symbols represent the observational LFs derived by dif- 
ferent authors up to z ~ 1.2 (upper panels of each figure). We 
have overplotted a PLE model in all the panels (dashed-dotted 
lines). The abscissas axis (absolute magnitude axis) in Fig.[3]has 
been transformed into stellar masses of m ETGs at z — using 
the expression derived for local ETGs by ICimatti et al.l {2006) 
and considering the modelled L-evolution of mETGs (axis at 
the top of the figure). The vertical dashed lines in panels a-e 
of Fig.|3]indicate the magnitudes that correspond to ETGs with 
log(A1*/A1o) = 11 and 12 at z = 0. So, our results concern to 
ETGs with M, > 10 n A1 o at z = 0, and roughly to L > V 
galaxies at z < 1 . 

Attending to the LFs of red galaxies at z < 0.8 (panels a to 
c of Fig. [3), we corroborate the negligible number evolution that 



all the observational studies report up to this redshift for the lu- 
minous red galaxies: notice that the PLE model for ETGs lies 
exactly on top of the observational data. However, observational 
studies disagree on the number evolution experienced by red 
galaxies at 0.8 < z < 1 (see panel d in the figure). Depending on 
the author, the number density of ETGs is found to decrease by a 
factor of ~ 2-3 or negligibly. The question is why authors derive 
such a different number evolution for red galaxies at 0.8 < z < 1 . 

The predictions of our model assuming that only ETGs are 
classified as red galaxies are also shown in Fig.[3](thin solid lines 
in top panels with lined area below them). The model is capable 
of reproducing the negligible number evolution reported by ob- 
servations for massive, red galaxies up to z ~ 0.8 (panels a-c in 
the figure). However, it predicts a significant drop in the num- 
ber density of mETGs at 0.8 < z < 1 (panel d), reproducing the 
data for red galaxies obtained by F07 at that redshift bin. So, the 
model shows that it is completely feasible to build up ~ 50-60% 
of the present-day number density of mETGs through the major 
mergers reported by observations since z ~ 1 . Another striking 
result of the model is that it predicts that the bulk of this recent 
mETGs assembly takes place basically during a time period of 
~ 1 Gyr elapsed at 0.8 < z < 1, in agreemen t with observa- 
tions by F07, but in contrast with other studies (IBell et al.l2 004b: 
ICool et al.tf 2008). However, the contradiction between these ob- 
servations and our model predictions is just apparent, as shown 
below. 

Let us now consider that DSFs are also identified as red ob- 
jects. In Fig. [3] we have also overplotted the LFs of ETGs+DSFs 
in the panels corresponding to the red galaxy population (thick, 
solid lines). Notice that, if DSFs are included into the red galaxy 
sample, the model now reproduces the data reporting an appar- 
ent neglig ible number evolution of massive, red galaxies up to 
z ~ 1 (Bell et al.l l2004b: Co ol et al.1 12008). This result suggests 
that the discrepancy on the number evolution reported by differ- 
ent studies for bright, red galaxies up to z ~ 1 could be due to 
the inclusion of a significant amount of DSFs into the red galaxy 
sample, depending on the color criteria used by each author. 

This possibilit y has also be e n suggested by other re - 
cent studi es (Bram mer et al.l 120091 ; iBenson & Devereuxl 12009). 
iKavirai et al.l (l2009al) reports that, accounting also for the late- 
type progenitors of mETGs in a ACDM scenario (and not just 
those of early-types), more than 50% of the stellar mass which 
ends up in mETGs today is in late-type progenitors at z ~ 1 (in 
agreement with our results). Our model corroborates their con- 
clusion on the relevance of accounting for the progenitor bias 
in observational studies, in order to derive accurate conclusions 
regarding the evolution of mETGs in the last ~ 8 Gyr. 

Note that the model predicts a contamination of the galaxy 
red sequence by dusty galaxies amounting to ~ 50% at 
z ~ 1 for our magnitude range, in agreement with observa- 
tions (ICimatti et al.l2002l;lYan & Thomp son 20031; IGilbank et al.1 
120031; iMoustakas et al.ll2004l) . This means that the existence of 
this DSF population at z ~ 1 can be explained just accounting 
for the gas-rich major mergers strictly reported by observations. 

The model starts failing to reproduce observational data at 
z > 1.2, independently on if DSFs are classified as red galax- 
ies or not. In 94.31 we discuss in detail the model limitations, 
especially at redshifts z > 1 . 

Our scenario also works for disk galaxies. Attending to the 
panels a-e of Fig.|4] we can see that the model predictions for 
blue galaxies match their observational evolution at z < 1 .2 as 
reported by F07, better if DSFs are considered as blue galaxies 
together with Sa-Irr's (thin, solid lines in the panels) than if not 
(thick, solid lines). The model predicts naturally the decrease 



9 



Eliche-Moral et al.: Reconciling the hierarchical buildup of mETGs at z < 1 with mass-downsizing 



log[M(ETG@ z = O)/M ] 

13 12 11 13 12 11 10 12 11 10 12 11 10 12 11 10 




-24-23-22-21-20 -24-23-22-21-20 -24-23-22-21-20 -24-23-22-21-20 -24-23-22-21-20 

M B (mag) 

Fig. 3. Evolution of the B-band LFs with redshift of red gala xies up to z ~ 1 .2, in redshift bins of Az - 0.2. Up per panels: 
Comparison with observational B-band LFs. Triangles: Data bv lBell et al.l d2004bl) . Crosses: Data bv lCool et alJ d2008l) . Pentagons: 
Data by F07. Lines: Model predictions. Thin, solid lines: Predictions considering that only ETGs are classified as red galaxies. The 
area below them in each panel has been lined in benefit of clarity. Dashed line: LFs at z — 0-0.2, according to the model. Dashed- 
dotted lines: Evolution of the LFs at each redshift in case of PLE. Dotted shaded region at the right of each panel: Magnitude 
range where the model is not valid, because galaxies have M(B) > My lm (B). Thick, solid lines: Predictions considering that also 
DSFs are classified as red galaxies. Lower panels: Model uncertainties, derived from the merger timescales (shaded area) and 
from errors in LSJ09 merger fractions (lined area). Typical error bars of observational data have been overplotted at three different 
magnitudes, covering the magnitude range of the model. Solid lines: Limits of the total uncertainty region. [A color version of this 
plot is available at the electronic edition]. 



by ~ 30-40% of the number density o f blue galaxies sin ce z ~ 
1 to z ~ reported by observations (Lib ert et al.\\200& F07), 
just considering the transformation of disks into ETGs driven 
by the major mergers at z < 1. Moreover, the majority of these 
gas-rich progenitors are predicted to be undergoing a transitory 
DSF phase of the major merger at 0.8 < z < 1 (compare the 
thick, solid line and the thin, solid line in panels i-j of the figure). 
Obviously, as the model reproduces the observational evolution 
of the LFs of blue and red galaxies up to z ~ 1, it also fits that of 
the total LFs (see panels a-e in F'ig . [5j> . 

In the lower panels corresponding to each redshift bin of 
Figs.[3][5] we show the model uncertainties due to the merger 
timescales and to the errors in the observational merger fractions 
(see jj2.5l ). These uncertainties are relevant enough to make diffi- 
cult to draw a conclusion on the production rate of ETGs, as they 
could imply from no number evolution of mETGs to a decre- 
ment of their number density by a factor of ~ 4. However, we 



must keep in mind that the nominal values used in the model are 
robust and realistic, derived as averages of the possible obser- 
vational values. By the contrary, the values used here for uncer- 
tainty estimates are quite extreme. Nevertheless, we will show 
that these uncertainties affect more weakly to the results in the 
/f-band ( 33.1.21 ). allowing us to rule out the PLE scenario not 
only for red galaxies at 0.8 < z < 1, but also for the blue galaxy 
population. This supports the robustness of the results derived 
for the B-band in this section. 

Summarizing, the model offers a framework in which the 
studies deriving a negligible evolution of red galaxies up to z ~ 1 
and those finding a significant evolution for them can be recon- 
ciled, just having into account that red galaxy samples can be 
highly contaminated by DSFs, probably intermediate stages of 
gas-rich major mergers evolving into ETGs. The model suggests 
that the follow up of populations selected by color (red vs. blue) 
backwards-in-time, as the majority of studies do since the pi- 
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Fig. 4. Evolution of the B-band LFs with redshift for blue galaxies up to z ~ 1 .2, in r edshift bins of Az = 0.2. Upper panels: 
Comparison with observational B-band LFs. Pentagons: Data by F07. Diamonds: Data bv lllbert et alj (120061) . Lower panels: Model 
uncertainties, derived from the merger timescales and from errors in LSJ09 merger fractions. The legend is analogous to that of 
Fig. |3] [A color version of this plot is available at the electronic edition]. 



oneering study by Strateva et al. (2001), might be providing a 
bi assed view of the gala xy assembly at z < 1, as already noted 
bv lKavirai et al.1 d2009ah . 

3.1 .2. K-band LFs for red and blue galaxies up to z ~ 1 

In Fig. [6] we show the model predictions on the evolution of the 
LFs for red, blue, and all galaxies for the Zf-band. The legend 
is the analogous to the one in Fig. [3] We have overplotted in the 
figure the model uncertainty region due to errors in the merger 
timescales and in the observational merger fractions by LSJ09. 

Comparing the observational data of red galaxies with the 
PLE model (panels g-i in the figure), we can see that red galaxies 
do not show any relevant number evolution at 0.25 < z < 0.75 
(panel g). However, their number density decreases suddenly by 
a factor of ~ 2-3 at 0.75 < z < 1, laying quite below the PLE 
model (panel h). So, ^T-band data are ruling out directly a PLE 
scenario for bright, red galaxies at 0.75 < z < 1. 

We have overplotted the Zf-band LFs predicted by the model 
in the same panels. Assuming that exclusively ETGs are iden- 
tified as red galaxies, the model predicts the negligible number- 
evolution observed for this galaxy population at 0.25 < z < 0.75, 
as well as the significant drop in their number density observed 
at 0.75 < z < 1 (see panels g and h in Fig.|6]l. The good agree- 
ment of the model predictions with Zf-band observational data 
even considering the model uncertainties stresses the feasibility 



of building up ~ 50-60% of the present-day number density of 
mETGs accounting for the major mergers reported by observa- 
tions at z < 1 • Again, the model predicts naturally the observa- 
tional fact that the bulk of this assembly occurs during the pe- 
riod of ~ 1 Gyr elapsed at 0.75 < z < 1, as in the Z?-band (see 
§13.1. 11 1. Moreover, considering that all DSFs are classified as red 
galaxies (thick, solid lines in panels g-i in Fig.|6|l, the model still 
reproduces better the observational results on red galaxies than 
the PLE model. 

Panels d-f in Fig.[6]represent the observed LF evolution for 
blue galaxies in the A'-band up to z ~ 1 .25. The figure shows that 
the PLE model (dotted-dashed lines) lies quite far away from the 
observational data not just for red galaxies, but also for the blue 
ones. So, the PLE scenario for blue galaxies is also discarded 
just considering observational /T-band data. 

The predictions of our nominal model on blue galaxies are 
overplotted in the panels. If DSFs are not included into the blue 
galaxy population (thick, solid line in panel d), the model for 
blue galaxies (Sa-Irr's) overlaps with the PLE model at 0.25 < 
z < 0.75, not fitting data. If the DSFs predicted by the model are 
considered as blue galaxies (thin, solid line with lined area below 
it), the model reproduces the data better than before, although it 
still underestimates the number density of blue galaxies at this 
redshift bin. However, DSFs become essential to reproduce the 
huge increase in the number density of blue galaxies reported by 
observations at 0.75 < z < 1, with respect to their present-day 
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Fig. 5. Evolution of the B-band LFs with redshift for all galaxies up to z ~ 1.2, in reds hift bins of A z — 0.2 . Upper panels: 
Comparison with observational B-band LFs. Pentagons: Data by F07. Squares: Data by IScarlata et alj d2007al) . Lower panels: 
Model uncertainties, derived from the merger timescales and from errors in LSJ09 merger fractions. The legend is analogous to that 
of Fig. |3] [A color version of this plot is available at the electronic edition]. 



value (by a factor of ~ 10, compare our model with the PLE 
model in panel e of Fig.[6]l. So, the model predicts two obser- 
vational excesses of blue galaxies with respect to PLE models 
simultaneously: 1 ) that of bright blue galaxies registered in the 
B-band, and 2) that of massive blue galaxies detected in the K- 
band. And both are explained just considering the gas-rich pro- 
genitors of the mETGs that have been assembled through major 
mergers at z < 1 . As in the B-band, the majority of these gas-rich 
progenitors are predicted to be undergoing a major merger at 
0.8 < z < 1. The model reproduces the evolution of total /if -band 
LFs up to z ~ 1 .25 better than a PLE model, even accounting for 
the model uncertainties (see panels a-c in the figure). 

3.1 .3. /-band LFs for red and blue galaxies up to z ~ 1 

In Fig. |7] we compare the LFs for red galaxies predicted by the 
model in the F814W-ba nd (~ /-band) with the observational 
results bv llm et aD dl996l) . The model reproduces quite properly 
the data at z < 0.8. However, notice that data at 0.8 < z < 1.2 are 
better reproduced if DSFs are being identified as red galaxies. 

3.2. LF evolution per morphological types upatz< 1 

In Fig.[TJ we plot the predicted evolution with redshift of the 
B-band LFs of different galaxy morphological types. DSFs 
have been included into the LFs of the Sc-Irr+M type, as 



they represent stages of major mergers with strong morpho- 
logical distortions. Observati onal data by different authors are 
also shown for compariso n dllbert et alJ2006l : lFranceschini et al.1 
2006; IZucca et al] |2006). The legend for lines is analogous to 
that of Figs. 061 but now refers to morphological types, not to 
color classes. 

FigureQ] shows that the model can reproduce the observed 
evolution of the LFs for different morphological types at z < 1 : 
E-SOa's (panels a-e), Sa-Sab's (panels f-k), Sb-Sbc's (panels 1- 
r), and Sc-Irr's (panels s-y). The figure suggests again that the 
major mergers that are observed at z < 1 provide a natural ex- 
planation for the observed numerical growth of mETGs since 
Z ~ 1, and that the gas-rich progenitors of these recently-formed 
mETGs can explain the observational excess by a factor of ~ 4-5 
of late-type galaxies at 0.8 < z < 1 with respect to a PLE model 
(see panel w in the figure). Intermediate types (Sa-Sab's and Sb- 
Sbc's) are consistent with a PLE evolution up to z ~ 1, according 
to their typical SFHs. 

Figure|4]shows that the decrease in the number density of the 
spiral (i.e., blue) galaxies from z ~ 1 to z ~ is ~ 30-40% of 
its present-day value. However, the same numerical decrease in 
terms of the present-day Sc-Irr+DSF population means that this 
population was higher by a factor of ~ 4-5 at z ~ 1 than now 
( Ilber t et al.1 l2006tl Francesch ini et all 120061: 



Zucca et alJ 2006; 



Bro wn et al.1 l2007t IScarlata et al.ll2007bt ICool et all 120081) . So 

the scenario proposed in the model for describing the evolution 
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Fig. 6. Evolution of the /T-band LFs with redshift of red galaxies (panels g-i), blue galaxies (panels d-f) and all g alaxies (panels 
a-c), a s predicted by the mo del, up to z ~ 1 .25. S ymbols: Observational LF data in the /T-band. Diamonds: Dat a bv [ Cirasuolo et al.l 
d2007l) . Hexagons: Data bv iPozzetti et al.1 (l2003h for 0.2 < z < 0.65 and 0.75 < z < 1.3. Asferis/b: Data bv lArnouts et alJ (120070 
for 0.8 < z < 1.0 and 1.0 < z < 1.2. Lines: Model predictions. The legend for lines is the same as in Fig. [3] but now for the 
/f-band. Shaded area: Total uncertainty region due to errors in the nominal values used in the model for the merger timescales and 
the merger fractions. Typical error bars of observational data have been drawn along the whole magnitude range in all the panels, 
just as reference. [A color version of this plot is available at the electronic edition]. 



of the galaxy LFs selected by color up to z ~ 1 also explains the 
observed evolution of the LFs selected by morphology. 

4. Discussion 

4. 1 . Mass-downsizing and the hierarchical assembly of 
mETGs at z < 1 

The most striking result of the model is that it reproduces the ob- 
servational evolution of the LFs at z < 1, predicting at the same 
time a late hierarchical assembly for at least half of the present- 
day mETGs (at z < 1). Recent semianalytical model within the 
standard ACDM framework reproduce the observed increase of 
the number density of mETGs since z ~ 1 (Ho pkins et alf 2008: 



iKavirai et al.l2009aUOesch et al.ll2 009: Catt aneo et al.l2010t) . As 
our model proves that the major mergers strictly reported by cur- 
rent observations are capable of explaining the mETGs assembly 
since z ~ 1 , our approach is directly linking the observed number 
of major mergers and the predicted one by hierarchical models, 
demonstrating that the hierarchical scenario and observations are 
quite compatible afterwards. 

Studies supporting galaxy mass-downsizing report that only 
~ 40-50% of the galax ies with At, > 10"A4 q were completely 
assem bled at z ~ 1 dPerez-Gonzalez et al.l 120081 lllbert et all 
2010). Having into account that ETGs dominate the mass func- 
tion at this mass range, we can conclude that the main result 
of our model is in agreement with the global trends of mass- 
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Fig. 7. Evolution of the /-ba nd LFs for red g alaxies with redshift, as predicted by the model. Pentagons: Observational /-band LFs 
for red galaxies derived bv llm et all (Il996h . Lines: Model predictions. The legend for lines is the same as in Fig. [3] but for the 
/-band. Shaded area: Total uncertainty region due to errors in the nominal values used in the model for the merger timescales and 
the merger fractions. [A color version of this plot is available at the electronic edition]. 



downsizing. Moreover, the model predicts that mETGs have ac- 
quired their actual volume density at z ~ 0.8, again in agree- 
ment with the previous studies. In Paper II of this series, we 
will show that the assembly epoch derived by the model for field 
mETGs is in agreement with that derived by semianaly tical mod- 
els and studies on th e SPs of mETGs (see De Lucia et al.l F2006; 
iNeistein et al.ll2006). Moreover, the model predicts that ~ 50% 
of the present-day number density of mETGs already exist at 
Z ~ 1, consistently with the detection of numerous "read&dead" 
mETGs in high-density environments at higher redshifts (see 
references in HJ. 

The model assumes that the numerical growth for each 
mass bin at z > 0.6 is similar, in agreement with several 
observations (see assumption #7 of 32.3.2b . This means that 
the mETGs at different masses are assembled at the same nu- 
merical rate at 0.6 < z < 1 in the model. Obviously, this 
growth rate is in contradiction with mass-downsizing, that de- 
rives that the most-massive galaxies must have experienced a 
lower number evolutio n at z < 1 than galaxies with lower masses 
dCirasuolo et all 12007b IScarlata et alll2007allbh . In fact, ETGs 
with mass es logfAL/Atp) > 11.5 seem to be already formed 
at z > 1 .7 (Frances chini et al1l2006h . This does not mean that the 
model and mass-downsizing can not be reconciled. Firstly, be- 
cause it can not be ruled out that mass-downsizing can be an arti- 
fact of large observational errors at t he bright-end of the de rived 
galaxy mass functions, as claimed by Fonta not et alJ (l2009). And 
secondly, because even if galaxies more massive than lO^Af© 
were already in place at z > 2, the model could reproduce it just 
considering a slower numerical growth rate at the high-mass end 
than at lower masses, without affecting its main results (these 
galaxies represent < 5% of the total mETGs population consid- 
ered by the model). However, this would require the prior knowl- 
edge of the distribution of major mergers with galaxy masses up 
to z ~ 1 , a question unsettled by the moment. 

In conclusion, the model shows that the observational phe- 
nomenon of mass-downsizing can be compatible with the hierar- 
chical assembly of nearly half of the present-day number density 
of mETGs at z < 1 . In fact, other recent studies confirm upsizin g 
in galaxies at z < 1.5 ( see iRvan etal]|2007l: fFranx et alj|2008l) . 
However, a complete reconciliation of hierarchical models and 



mass-downsizing might require a better implementation of the 
physical processes affecting to massive systems in the simula- 
tions, and a more exhaustive analysis of the biasses and uncer- 
tainties of observational data. 

4.2. Comparison with other observational results 

The model offers a framework for galaxy evolution that can 
explain or agrees with several recent observational results: 

1. The bright-end of the LF in Ha at z ~ 0.84 is dominated 
by Ir r and merging systems (Willar etal. .1 \2008[ \Sobral et al.\ 
2009). — The model predictions are coherent with this fact, as 
the gas-rich progenitors at 0.8 < z < 1 of the recently-assembled 
mETGs represent ~ 30-40% of the present-day number density 
of all spirals, and are Sc-Irr systems. Therefore, the model 
predicts that ~ 50% of the star-forming galaxies at z ~ 0.8 
are Sc-Irr or DSF merging systems. As Sc-Irr's have enhanced 
Ha emissions as compared to normal disks (and much more i f 
they are involved in a major merger, see West raetalJ 12010). 
the model predicts that they must control the Ha luminosity at 
z ~ 0.8, as observed. 

2. Nearly 75-85% of quiescent galaxies with 
logfAL/ Atp) > 10 11 at z < 0.8 have an elliptical mor- 
phology filbe rt et al.\ \2010\) . — The model predicts that the 
major mergers start affecting the LF of ETG appreciably at 
z > 0.8. At lower redshifts, the mETGs are nearly in place in the 
model (some evolution is still ongoing, but it is negligible, see 
figures in §|5}, as observed. 

3. An increased fraction o f the SF takes place in ETGs at 
z > 0.8 ftDahlen et al.\ \200^) . — The model predicts that the 
bulk of the number evolution of ETGs through major mergers 
at z < 1 basically takes place during a period of ~ 1 Gyr at 
0.8 < z < 1. As these major mergers are predominantly wet, 
their remnant mETGs must have undergone strong starbursts 
during the merging-nuclei phase of the encounter, meaning 
that a relevant fraction of the SF at z > 0.8 must be occurring 
in spheroidal remnants that are evolving into typical ETGs, in 
agreement with the previous result. 



14 



Eliche-Moral et al.: Reconciling the hierarchical buildup of mETGs at z < 1 with mass-downsizing 



4. Redshift z ~ 0.8-7 is a transition epoch in the formation 
history of ETGs ACimatti et al.\\2d06\ \Franceschini et al . 2006; 
\llbert et al]\201(A) . — The number evolution of mETGs is nearly 
frozen at z < 0.8 in th e model, in agreement with numerous 
observational studies ( van der Wei et al. 2007 ;_ Banerii et alj 



120101 iHuertas-Companv et alJ l2009a: Rudni ck et al. 2009). The 
slow down of the generation of ETGs at z < 0.8 is due to two 
facts: basically, to the observational decrease of the merger 
fractions at these redshifts and, secondarily, to the more relevant 
role played by dry and mixed mergers at z < 0.8 in comparison 
to wet events (see Fig.|2]i. Our model proves that the declining 
of the wet mergers rate at z < 1 can explain the observed slow 
down in the assembly of mass ive ellipticals since that epoch, as 
proposed bv lllbert et alj d2010l) . For more details, consult Paper 

n. 

5. The evolution of the global LF is mainly driven by the SF 
associated t o disk galaxies, inste ad of by the number evolution 
of galaxies AWalcher etal.\\200$) . — If we compare the LFs of 
the total galaxy populations with the LFs of disks and ETGs in 
Figs.[3]|5] we can conclude that the evolution of disk galaxies 
(which is basically due to their SFH) determines the global 
shape and behaviour of the LF of all gal axies at z < 1, in 
agreement with this result bv lWalcher et alJ 

6. The SFRs of star-forming galaxies seem to be largely 
unaffected by the loc al processes that truncate SF at z ~ 0.8 in 
high density regions iPatel et ~al. 2009) . — In the model, the SF 
quenching takes place during the post-merger phase of major 
mergers. Normal disks not involved in major mergers continue 
with their exponentially-decaying SFR, without experiencing 
SF quenching. If we also consider that the enhancement of 
the SFR due to major mergers is negligible as compared to 
the typical SFR un d ergone by normal di s ks even at z ~ 0.8-1 
(Sobr aletal.1 l2009t iRobaina et all l2009t lLopez-Saniuan et al l 
2010), the model implies that the global SFR of typical disks is 
not affected by the processes that shut down the SF at z ~ 0.8, 
in agreement with the result stated above. 

7. Observational major merger rates must account for 
~ 50% of the mass build-up at the massive end at z < 1 
iWalcher et al.\ \200&) . — At the brightest magnitudes, where 
the LF of mETGs dominate over the rest of galaxy types in 
the global LF, the number evolution of mETGs determines the 
evolution of the global LF, and explains the buildup of ~ 50% 
of the systems in this range (see Figs.[3]|6]l. The late buildup 
of mETGs is being confirmed by recent spectral studies that 
conclude that ~ 50% of the total stellar mass density of mETGs 
has appeared at z < 1 (iThomas et al.l 12005b iFra nceschini et al.1 
2006- ICool etaLll2008Hllbert et al.ll2010t). Our model proves the 
feasibility of the conclusion derived bv lWalcher et al.l 

8. A signif i cant DSF population appears a t 0.8 < z < 1 

( Cimatti et al\ ~2002\ Yan & Thompson] \2003t \Gilbank et al.\ 
120031 - 1 Moustakas et al. 2004; Weiner et al. 2005). — The con- 
tribution of gas-rich major mergers is negligible at z < 0.7, and 
thus DSFs (understood as transitory stages of major mergers 
by the model) are not relevant n umerically at thes e redshifts, 
in agreement with observations (Bell et all l2004al) . However, 
the model predicts that large amounts of DSFs must appear 
at z > 0.8, associated to the rise of gas-rich major mergers. 
According to the model predictions, ~ 50% of the galaxies 
in the massive-end of the red sequence at z ~ 1 must be 
DSFs, in agreement w i th observations dCimatti et aT 20021 



9. The build up of ETGs at z < 1.2 seems to have been 
confined only to low-density en vironments, and not to galaxy 
clusters \Pozzianti et al.\\2009b\) . — The model is also coherent 
with this observational fact. Notice that it does not consider 
the effects of the environmental density on galaxy evolution, 
but a global, numerical evolution of field galaxies. In fact, the 
observational merger fractions and local LFs that we use are 
derived from samples of field galaxies. So, the model is tracing 
the average build up of mETGs in the field since z ~ 1 . 

10. Above AL > 1O 1O - 8 A1 , blue E-SO's at z < 1.4 resemble 
to merger remnants pro b ably migrating to the red-sequence 
\Huertas-Companv et al\ 1207 Ol) . — Our model is coherent 
with this result, and proves the feasibility of reproducing the 
observed mass migration from the blue galaxy cloud to the red 
sequence atz < 1 for masses At, > 10 1I A / l o , just accounting for 
the major mergers reported by observations. 

1 1 . The evolution of luminous red galaxies since z ~ 
0.5 is consistent with t hat expected from passive evolution 
I Toieiro &■ PercivaA\20lb]) .— Our model is in agreement with 
this result, as the buildup of mETGs is nearly frozen at z < 0.7 
(see paper II of this series for more details). 

4.3. Model limitations 

The model reproduces the evolution of the number density of all 
galaxy types at z < 1, but fails at z > 1 (see Figs.[3][T]). Several 
aspects on galaxy evolution at z > 1 can be limiting the validity 
of the model at higher redshifts: 

i. Validity of the Hubble classification scheme at z > 1 ■ — 
It is known that the fraction of galaxies with highly-distorted 
morphologies increases at z > 1 . In fact, the e mergence of the 
Hubble Sequence seems to occur around z ~ 1 .4 (lOdew ahn et al 
1996; IWu et alJl200ll: iPapovich et all 120051: iRavindranath et al 
2006). This means that the correspondence between SFH and 
morphological type (as assumed in the model) lose its meaning 
atz > 1. 

ii. Significant contribution of other evolutionary 
processes. — Besides major merging, other processes are 
capable of transforming spirals into SO-SOa galaxies, as the 
effects of the intracluster medium on the infallen star-forming 
galaxies or the minor mergers (Eliche-Moral et al. 2006a 



| Aragon- S alamanca 120081 
2009a UGeach et alJl2009t 



Cassata et al. 2008; Poggianti et al. 



Temi et all 1200 9: Wolfetal. 2009 



However, it seems that the effects of these processes are 
more relevant in the formation of intermediate-mass disks 
of types later than Sa than in the buildup of ETGs at z ~ 1 
dAzzollini et alJl2009HWeinzirl et al.ll2009tlL6pez-Saniuan et all 
2010). At higher redshifts, the relative relevance of these 
processes might be quite different. In fact, observations at high 
redshift and cosmological simulations favour a scenario in 
which galaxies at z >> 1 formed through steady, narrow, cold 
gas streams that pe netrate the shock-heated media of massi ve 
dark matter haloes dlacconi et alj 120081 iDekel et all [2009aHbl: 
Fors ter Schreiber et alJl2009tlCeverino et al 

Hi. Higher gas amounts in disks at z > 1. — Major mergers 
at z > 1 basically involve gas-rich gal axies with higher gas 



Yan & Thompson! 120031 iGilbank et al.l 120031 iMoustakas et all 
2004l lWild et alJl2009blah . 



contents than thei r local counterparts jPapovich et al. 2005 
| Genzel et al.l 120081; iTacconi et al.l 120081 lF5rster Schreiber et al. 
2009; lLaw et al.l 120091) . Recent observations and simulations 
show that major mergers in extremely gas-rich media can 
give place to galaxies of later types tha n SO, as gas infall 
can rebuild the disk in ~ 2-3 Gyr (see ISchweizer & Seitzeri 
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2007t iGovernato etaDl2009l: lHammer et af] l2009b). So, as this 
situation seems to be common for major mergers at z > 1, 
the main model hypothesis may be not fulfilled, invalidating 
model results at z > 1 . Moreover, numerical simulations suggest 
that extremely gas-ri c h disks tend to form massive clumps 
ddi Matteo et alj|2008t lElmegreenll2009t) . These clumps could 
lead to an overestimation of the real fraction of major mergers at 
z > 1 through the classical CAS techniques, making the model 
predictions unrealistic. Higher gas reservoirs in disks at z > 1 
can also imply more massive starbursts induced by mergers. 
In this regard, CAS techniques could identify a minor merger 
event as a major one, overestimating again the number of major 
mergers at z > 1 . 

iv. Observational merger fractions at z > 1. — The merger 
fractions used in the model are robust at z < 1 , where they were 
derived (see 92.3.21 1. Results for higher redshifts use the extrap- 
olation of the trend at z < 1 to earlier epochs. Although this 
extrapolation is very similar at 1 < z < 1 .5 to the merger frac- 
tions derived bv lConselice et al] (20081). .they may be overesti- 
mated by a factor of ~ 2-3 djogee et aDl2009h . contributing to 
the uncertainties of the model results at z > 1 . 



5. Summary and conclusions 

We present a model that predicts the backwards-in-time evolu- 
tion of the local mETG population (with \og(M„/M Q ) > 11 at 
z = 0), under the hypothesis that each observed major merger 
gives place to an ETG. The model traces back-in-time the evo- 
lution of the local galaxy populations considering the number 
evolution derived from observational merger fractions and the 
L-evolution of each galaxy type due to typical SFHs. The main 
novelty of the model is the realistic treatment of the effects of 
major mergers on the LFs. In particular, we have had into ac- 
count the different timescales of a major merger, the colors of 
merging galaxies at each merger phase, and the different mor- 
phological types of the progenitors depending on whether the 
merger is wet, mixed, or dry. All the model parameters and evo- 
lutionary processes considered in the model are based on robust 
results of observational and computational studies. 

The model is capable of reproducing the observational evolu- 
tion of the bright-end of the galaxy LFs at z < 1, simultaneously 
for different rest-frame bands (B, I, K) and selection criteria (on 
color, on morphology). Moreover, it can explain several recent 
observational results, such as the morphology of the star-forming 
and quiescent populations at z ~ 0.8, the global evolution of the 
SFR at z < 1, and the fraction of red sequence contamination by 
DSFsuptoz~l. 

The model predicts the appearance of ~ 50-60% of the 
present-day number density of mETGs since z ~ 1 (in agreement 
with some observational studies), just accounting for the major 
mergers strictly reported by observations during this epoch. The 
bulk of this evolution has taken place during a time period of 
~ 1 Gyr elapsed at 0.8 < z < 1, also in agreement with obser- 
vations. So, the model reproduces the progressive migration of 
stellar mass from the bright-end of the blue galaxy cloud to that 
of the red sequence observed during the last ~ 8 Gyr through 
mixed and wet major mergers, and from the lower masses to the 
higher masses of the red sequence through dry and mixed merg- 
ers. In this sense, the model corroborates the feasibility of the 
mixed evolutionary scenario proposed by F07, al least for the 
bright-end of the LFs at z < 1 . 

The framework proposed by the model can also reconcile 
conflicting observational results on the amount of number evo- 



lution experienced by massive red galaxies, just considering the 
usual contamination by DSFs exhibited by red galaxy samples 
at z > 0.7. Moreover, the number evolution of massive galax- 
ies predicted by the model at z < 1 is compatible with global 
trends of mass-downsizing, a fact that suggests that the hierar- 
chical scenario and this observational phenomenon can be made 
compatible. 

The model provides several secondary results on the evolu- 
tion of galaxy populations, summarized below: 

- The model reproduces naturally the excess of bright late-type 
galaxies (by a factor ~ 4-5) observed at 0.8 < z < 1 with re- 
spect to PLE models (or, equivalently, of ~ 30-40% of bright 
blue galaxies), just considering the gas-rich progenitors of 
these recently-assembled mETGs. 

- It also explains the observational fact that z ~ 0.8 is a tran- 
sition redshift in the recent assembly history of mETGs. The 
combination of the decrease of major merger fractions at 
z < 1 with the dominant role played by mixed and dry 
mergers since z ~ 0.8 over wet mergers makes the assem- 
bly of mETGs to be nearly frozen at z < 0.8 (in terms of net 
growth). 

- The model predicts the appearance of a significant popula- 
tion of DSFs at z > 0.8, associated to the rise of gas-rich ma- 
jor mergers at these epochs. The predicted number density 
of DSFs at z ~ 1 (considered as transitory stages of gas-rich 
major mergers) can explain the ~ 50% of contamination of 
the red sequence by dusty galaxies reported by observations 
at z ~ 1 • This DSF population is essential for reproducing 
the noticeable observational number evolution experienced 
by bright blue galaxies in the /T-band since z ~ 1 . 

The role of major mergers in the assembly of mETGs has 
been supported by several observational studies, but this is the 
first time that it is proved the feasibility of building up nearly 
half of the present-day number density of mETGs through the 
major mergers strictly reported by observations at z < 1. The 
model provides an evolutionary scenario that establishes natu- 
rally a link between the major mergers occurring at z < 1, the 
appearance of a significant population of DSFs at 0.8 < z < 1, 
and a significant hierarchical assembly of mETGs during the last 
~ 8 Gyr, in agreement with predictions of hierarchical models 
and with observations. 
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Appendix A: Decomposition of the LFs in r* by N03 
in finer morphological type bins 

A.1. General description of the procedure 

As commented in 92.2.2I we have decomposed the r*-band LFs 
derived by N03 into thinner classes, using the information of the 
LFs in the B-band by C03. The morphological types defined by 
C03 are thin enough as to reproduce each of the coarser mor- 
phological classes defined by N03 as a sum of "thin" types. 
As the morphological classifications of a galaxy in the B and 
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r* bands are similar for local galaxies (Kuchi nski et al.l 12000). 
we can expect that a galaxy at z ~ exhibits nearly the same 
morphological classification at different bands. If a univocal re- 
lation between the absolute magnitudes in the B and r* bands 
exists for each galaxy type at z ~ 0, the number of galaxies of 
a certain type with magnitude M(r*) would be the same as the 
number of galaxies of the same type exhibiting the correspond- 
ing M(B) magnitude. Let us assume that this univocal relation 
exists at z ~ 0, at least fulfilled by all galaxies within the same 
morphological type (we will justify it in {3A.2I ). This means that, 
for all the galaxies within a given morphological type: 

M(B) = flo + a\M(r*), (A.l) 

being Oq and a i constants that depend on the morphological type. 
We can also assume that all the galaxies detected in B are going 
to be also detected in r* (this is true at least for bright magni- 
tudes). Therefore, the number of galaxies with M(r*) and with 
M eq (B) will be the same for each morphological typ4B In partic- 
ular, this occurs for the E's and SO's, and hence: (p[M(r*), "E"] = 
0[Mg q (B), "E"] and <p[M(r*), "SO"] = 0[M^,(B), "SO"]. Notice 
that the M(B) magnitude corresponding to the same M(r*) for 
different morphological types will be different a priori (that is 
why we note them by the subindices "E" and "SO"). The previ- 
ous expressions imply that: 

<p[M(r*),"E"]/<l>[M(r'),"S0"] = 0[M= q (5),"E"] /0[M= q o (5),"SO"].(A.2) 

This equation would be valid for all M(r*). The last term of 
the previous equation is given by the LFs by C03, provided 
we know the M(r*)-M(B) relation for each morphological type. 
Obviously, we can consider that: 

(f>[M(r*\ "E - SO"] = <f>[M(r*), "E"] + (f>[M(r*), "SO"], (A.3) 

where the left term is given by the LF derived by N03. Therefore, 
if the univocal relation M(r*)-M eq (B) for each galaxy type ex- 
ists and it is known, we can decompose the local LF in r* for 
E-SO's into thinner classes (0[M(r*), "E"] and (f>[M{r*), "SO"]), 
using eqs. IA.2l and lA.3l 

As this procedure can be performed also with the rest of mor- 
phological classes, we could split the four LFs in the r*-band 
by N03 into the eleven thin types defined by C03, as required. 
Nevertheless, the introduction of so many types in the model 
does not insert any additional advantage, because the character- 
istic properties of adjacent, thin morphological classes are quite 
similar. However, C03 also derived the local B-band LFs for six 
coarser morphological classes (E, SO-SOa, Sa-Sab, Sb-Sbc, Sc- 
Scd, and Sd-Irr). Note that the four classes by N03 are not strictly 
compatible with these six coarser types, in the sense that they can 
not be reproduced as simple sums of these six types, except for 
certain classes. Therefore, we have combined the information by 
N03 with the information for thinner and coarser types by C03 
to derive the LFs in the r*-band for the following morphological 
types: E-SOa, Sa-Sab, Sb-Sbc, Sc-Scd, and Sd-Irr. Although we 
are just considering only one morphological type more than N03, 
the new distribution of morphological classes groups galaxies 
with more similar spectrophotometric and morphological prop- 
erties than the classes defined by N03. The LFs in these five 
classes can be easily obtained through a procedure analogous to 
that used to derive eqs. IA.2l and lA.3l combining the LFs by N03 
and those by C03 in coarse and thin types. 

1 We will note the M(B) magnitude corresponding to a given value 
of M(r') through eq.|Aj]as M e \B). 



The resulting LFs are plotted Fig.Q] Their Schechter 
parametrizations are listed in TableQ] In order to obtain an ade- 
quate Schechter fit to the LF of ETGs obtained through this pro- 
cedure, we have fixed the a slope of ETGs. This does not affect 
the results in the model, as we are limited to the bright-end of 
all LFs by the limiting absolute magnitude of the observational 
merger fractions by LSJ09. The residuals between the total LF 
derived by N03 and the resulting one from the decomposition are 
lower than ~ 0.1 in M(r*) > -22 mag, peaking at the brightest 
magnitudes. Residuals between the LFs obtained for each mor- 
phological type are lower than 5% in the whole magnitude range 
of interest, supporting the robustness of the procedure. 

A.2. Proportionality between M(B) and M(r*) for local 
galaxies 

Zhao & Newber g| (120061) indicate that magnitudes in the B- 
band can be transformed to the r*-band, according to: r* = 
V - 0.451 (B - V) + 0.082. We have sele cted a local galaxy 
sample (z < 0.1) from the GOYA survey (lAbreu et al.l 12007). 
which basically contains spiral galaxies. We have checked that 
the (B - V) color distribution of local galaxies is quite narrow 
(< B - V >= 0.8529, with cr = 0.1218). So, for local spirals, we 
can deriv e that: r* ~ V - 0.303, by replacing this median color 
value into lZhao & Newbergl (12006) conversion equation. 

We have plotted the magnitude of these local galaxies in the 
V-band against their magnitudes in the B-band. A well-defined 
linear trend is found, with a Pearson linear correlation coefficient 
of p ~ 0.96. Therefore, we can consider that for local galaxies: 
B = c + c{V, being c = -1.3 ± 0.4 and c x = 1.02 + 0.06 valid 
for Sa-Irr galaxies. Replacing this relation into the previous r*-V 
relation, we find that: B ~ c\r * +(0.303ci + cq). This implies 
that, for local spiral galaxies, we can assume that m(B) oc m(r*). 
Considering that the k- and e-corrections at z < 0.1 are negli- 
gible and that global effects of dust-extinction can be neglected, 
the previous relation could be extrapolated to absolute magni- 
tudes, and thus: M(B) oc M{r*). Therefore, the assumption in 
eq. lA. 1 l is plausible for local disk galaxies, resulting for them: 

M(B) ~ 0.98M(r*) + 1 .6. (A.4) 

The existence of the narrow red sequence of galaxies in the 
local Universe implies that a similar relation must exi st for E- 
SO's, as it is basically composed by E-SO's at z ~ (Blanton 
l2006h . As the k- and e-corrections and the dust extinction are 
also negligible in the local Universe for E-SO's, we can consider 
that M(B) - M(r*) ~ c 2 also for E-SOa's. In order to obtain a, 
we have compared the LFs of ETGs in the B and r* bands, esti- 
mating the displacement on the magnitude axis required in order 
to make them to overlap. For local E-SOa's, we can derive that: 

M(B) ~ M(r*) + 1.5. (A.5) 

Equation lA.4l is consistent with the typical, observational colors 
derived for Sa-Sc galaxies (< B - R >~ 1.4), while eq. IA.5l is 
corrobora ted by typical colors of real E-S0 galaxies (< B -R >~ 
1.6, see lMannucci et aLl uOOl). Equations lA.4l and lA.5l prove the 
assumption adopted in eq. IA.ll The residuals obtained by com- 
paring the total LF derived by N03 in r* with the total LF by 
C03 transformed to the r*-band (using eq. lA.4llA~5l l are lower 
than ~ 0.01, supporting the assumption adopted in eq. IA.ll 
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Fig. 1. Predicted evolution of the B-band LFs by morphological types (E-SOa, Sa-Sab, Sb-Sbc, Sc-Irr) with redshift, up to z ~ 1.5. Observational d ata have been overpl otted for 
comparison. The abscissas axis (absolute magnitude axis) has been transformed into stellar masses using the expression derived for local ETGs bv lCimatti et al.l (120061) (axis at 
the top of the figure). Panels a-e: Evolution of the LF for E-SO's. Panels f-k: For Sa-Sab's. Panels l-r: For Sb-Sbc's. Panels s-y: For Sc-Irr+M's. Symbols: Observ ational data in 
the ij-b and for each morpho logical type. Typic al observational error bars in the whole magnitude range have been overplotted for comparison. Squares: Data by Scarla ta et al.1 
(2007a). Triangles: Data by IZucca et al.l d2006l) . Pentagons: Data for red galaxies by F07. Lines: Model predictions. The legend is analogous to that of Fig. [3] but referred to 
morphological classes, instead of color classes. Shaded area: Total uncertainty region due to errors in the nominal values used in the model for the merger timescales and the 
merger fractions. Notice that Sa-Sbc galaxies are not affected by major mergers in the model (and hence, they do not exhibit significant uncertainties). The uncertainties of the 
model predictions for ETGs are shown in Fig. [3] [A color version of this plot is available at the electronic edition]. 



